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A   STUDY  OF  THE  STRESS   DISTRIBUTION 
AROUND  ELLIPTICAL  HOLES   IN  A  FLATE  PLATE 
SUBJECTED  TO  IN-PLANE  LOADS 

by 
Demosthenes  HOZOS 


ABSTRACT 

A  theoretical  solution  of  the  title  problem  was  obtained 
and  a  computer  program  was  developed  to  examine  various  geometric 
arrangements  of  elliptical  holes  in  flat  plates  with  different 
loading  conditions. 

Graphs  of  theoretical  predictions  are  presented,  discussed 
and  compared  with  experimental  results  of  other  authors. 
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NOMENCLATURE 

2  C  Distance  between  the  foci  of  the  ellipses 

a ,  b         Major  and  minor  semiaxes  of  the  elliptical  hole 

I",  tj  Elliptical  coordinates  of  a  point  P(x,y) 

§"o  Constant  value  of  £"  of  any  point  on  the  elliptical 

hole 

^.Xl  '      Stress  functions 

<£*)  Normal  stress  intensity  in  the  tangential 

direction  to  §f  =.  constant 

65  Normal  stress  intensity  in  the  normal  direction 

to  %=   constant 

Tttf  Shear  stress  intensity 

The  usual  notation  for  stress  intensities  in 
Cartesian  coordinates 

Dimension  of  the  plate  in  the  x  direction 

Dimension  of  the  plate  in  the  y  direction 


6* ,  6j) 

i,  XM 

Wx 

Wy 

}x  = 

2CL 

Wx 

1y  = 

£b 

Wy 

Uniformally  distributed  load  applied  on  the 

edge  x=  — -  and  parallel  to  X  axis 
6a.  Uniformally  distributed  load  applied  on  the 

ed^c  y  =  —^  and  parallel  to  V  axis 
O-  ,  b        Semiaxes  of  the  confocal  ellipse  | -  |o  passing 

through  a  point  "PCx;y)^P(£>^') 
5  C  F        Stress  concentration  factor  =.  ■^rn<x:<.-  or 


SCF        Stress  concentration  facbor  =  5CF*  Cl-^)   or 
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u  =  bfa.  when  loading  parallel  bo  the  minor  axis  of  the 

'  elliptical  hole. 

Ll  -  cl/6  when  loading  parallel  to  the  major  axis  of  the 

elliptical  hole. 
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INTRODUCTION 

The  behavior  of  thin  flat  elastic  plates  perforated  with 
circular  holes  is  well  known  for  in-plane  loading. 

Theoretical  methods  and  experimental  results  have  been 
published  by  many  authors.  Refs.  (1-7) 

The  case  of  one  elliptical  hole  in  an  infinitely  large 
plate  is  discussed  in  reference  [l]  for  several  loading  conditions 
including  the  case  where  the  major  axis  of  the  ellipse  is  inclined 
to  the  axis  of  the  axial  load. 

Durelli  and  Murray  [10]  performed  photcelastic  tests  on  thin 
flat  plates  loaded  with  various  combinations  of  biaxial  stresses. 

Durelli,  Parks  and  Feng  [9]  examined  cases  of  rectangular 
plates  with  one  elliptical  hole  under  axial  loading,  again  using 
the  phot oe la stic  method. 

The  purpose  of  this  study  is  to  approach  the  problem 
theoretically  and  obtain  results  for  the  case  of  finite  thin 
rectangular  plates  perforated  with  one  elliptical  hole  and 
subjected  to  uniaxial  or  biaxial  loading. 

A  computer  program  was  developed  and  the  published  results 
of  available  experiments  were  used  for  comparison. 
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THEORY 

The  use  of  elliptical  coordinates  simplifies  the  mathematical 
analysis  of  the  stresses  around  an  elliptical  hole  in  a  flat  plate 
loaded  in  the  plane  of  the  plate. 

The  equations  of  transformation  from  Cartesian  to  elliptical 
coordinates  are 


X  =   C  cosk£  cos<q  (l) 


y  s  csiuh-l   sm^  (2) 


The  conditions  §  =  constant  and  ti  -   constant  are  represented 
"by  a  family  of  confocal  ellipses  and  a  family  of  confocal  hyperbolae 
respectively,  that  are  orthogonal  and  may  be  expressed  as 


X* 

+ 

ya 

=  i 

c^cosh1! 

X* 

ca  smK'f 

y2 

ca  cosx^i 

C2'  SW^V) 

0) 


00 


Figure  La.     represents  the  elliptical  coordinates  £",77 
of  a  point  P  (  x ,  y  )  . 
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Using  the  notation 


f-   _£  (co.sh2it-co5.2r?;)  (5) 


the  stress  components  can  "be  expressed  in  elliptical  coordinates  in 
terms  of  a  stress  function  Cp    as 


g         I      c)4>       sinhs"  cosh£"         c)  4> 


sinT7  CQST7  c)  4>  ^g) 


6? -  _L     ^^         sin77  cost?  c)4>     + 

J1     S72      Dx(s'm27j/  +  siTihY)      cJy 

SinhS"  cosh  £"  ^  4>  frj\ 


sinh^  cosh  £  o><^  ,3% 


Ja(srn*'n  +  si'nhl£)       ci-n 
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It  is  easy  for  one  to  verify  that  the  expressions  on  page 
538  in  [kl   reduce  to  equations  (6),  (7)  and  (8)  above. 

It  might  be  noted  here  that  the  following  assumptions  have 
been  made  in  the  foregoing  equations: 

i. Isotropic  elastic  material 

ii.  Thin  flat  plate 

iii.  Infinitesimal  deformations 

The  stress  functions  (^a-j)  assumed  here  are  the  same  as 
those  -which  Coker  and  Filon  used  in  reference  [h] 


Xl  -  e        cos  [0-0-7]  +  e         cos|>ui)«7j  (9a) 


X*  =  e        cosTo-iy']]  +e         Cos|Cn40'77j 
X5=  e^si^o-O^]  +  e(n"01  sm["(^07] 


(9b) 


(9c) 


X«  -  eCr,+°   Si*.[c*-0*Z]    +  e^"'     siu[(^H)^]  (9d) 

Xs  -  e       cos  ^7  (9e) 


Xg=    e.  C05^  (9f) 


X7  -  e1*  si-noi^  (9g) 
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Xa  =  e   SWi  ^  (9h) 


Xq=  5  (9i) 


X.o  =  «J  (9d) 


All  the  above  functions  satisfy  the  biharmonic  equation  for 
interger  (or  zero)  values  of  n.  I"  It  1 

Equations  (9a)  through  (9j)  can  be  combined  to  give  a  general 
expression  of  the  stress  function 


}a»  e    cos  (n-i)n  +  e    cos  (n  +  QT|   + 

r  fn+0*  .        cn-ol"  .       , 
4.  eni  e   sin(n-o^fe   sm(n  +  i)y  t 

r  r  -Cn+O?  •     ,    -(n-l)$  . 
-j-  fn|_€     Sin(n-i)i/4-e     Sin  (n-H)l?J  4. 

+   on   e    cos  ni^  4.  hn  e     cos  1177  4. 

4.  ]n    e     smrn^    4-ftn  e     Sin  r?-n  (10) 
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When  examinig  the  series  with  n=0  and  n=l  then  for 

finite  stresses  and  displacements  it  is  evident  that  the 

coefficients  d  ,f  ,h  ,j  ,k  ,f,  must  equal  zero. 
o,o,o,Jo'o,l       H 

The  stress  function  given  by  equation  (10)  can  be 
substituted  in  equations  (6), (7), and  (8)  for  the  stresses. 

However  in  the  present  analysis , only  loading  conditions 
that  maintain  symmetry  about  both  the  x  and  y  axes  are 
considered  and  therefore  the  coefficients  of  the  terms  that 
produce  nonsymmetric  functions  of  y   for  different  values  of  n 
were  set  eaual  to  zero. 

Finally  the  stresses  from  (6), (7)  and  (8)  have  the  form 


COS  (h+5>7  + 


+ 


Pn  (n+3)  € 


(n-3)f 


COS  (n4i)T|  4- 


+ 


^  Cne    -  Cn  (n-5)  e    -f  4  t>n  e 


-"Dn  (n-3)  e 


-(n+3)l" 


Cos  (n-l)^  4 


+ 


cmOl" 


Cn  (n-l)e    4  Dn  (n-\)  € 


-fn  +  0^  -i 


COS 


O7-O7  1  4 
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CO  >>■  v 

+  )_  n       &n(n-0€       +Hn(n-i)e      C05(n+2^7j?  4 


na2,4- 


+  [-&n(n-l)e         -Gn(n+Oe  -Hn(n-l)e 

-  Hn  (n  +  i)e  J   od5  nrj    4 

+  [Gn(n+Oe    4Hn(n+oe     JcosCn-27y     4 


4  2Ao  smh  2?" 


(11) 


00 


Cn-0$ 


-Cn-0$ 


,4  ^ —         r  -  Cn-OS  -<n-i;j   - 

C*  Tl  =  l,3...  (. 


n  +  i)? 


4    [-4Cn  €  4  Cn  (n-i) 


(n-3)?  -Cn4  0T 


4t)n(n-i)e  j  cos (n 41)17   4 


4[4Cne  +Cn(n+i")e  +  4Dn£  4 


4Dn  (n-fi)e  ]cosCn-r)7p  + 


n4i)| 


-Oi  +  iA 


[-Cn  (n  +  3)en         -Dn(n43)e  Jco5(n-3)^      + 


OO 


+  n3 


.n^ 


n|[-Gn(n-i)e     -Hn(n-0e    Jcos  (n+z)^    4 
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(n4;0f 


cn-r;s" 


+  [Gn(n-Oe        4  Gn(n  +  i)e         4Hn(n-i)e        4 
+   Hn  (ri4i)e  cosn^  4 

-  2  Ao  s  i  n.  h  2$ 


(12) 


4XV*    - 


:n-i)^ 


-Cn-Ot 


n[  [-CnCn-Oe""'  +'0n(n-0c  ]si*r?(n4-3)r>f 

4     [Cn(n4i)e  -T)nCn+0C  Jsin(in+0^4 


r  Cn43)£  -Cn+O*-,     . 

4    [Cn  (n-Oe         -PnCn-Oe  JsmCn-0^4 


Cn+ O^i 


r  Cn-fi)?  x-(n+i;S-,  ,      1. 


r  r  n^  ~n^  i 

t  11      [-GriCn-Oe   +Hn(n-i)e         sin(n+2)rj)  4 


4- 


Ga(n-0e  4-Gn(ri4i)e       -  Hn  (n-l)e 


-Cn-i)5 


-Hn  (n-n)e 


sin. 


n^  + 


-nr 


+  Hn  (n+0  e     Jsm  (n-ily  4  <?^o  Sin  ^r> 


(13) 
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However  equations  (11),  (12)  and  (13)  are  rearranged,  in  Appendix  1, 
in  a  more  suitable  form  for  translation  into  computer  language. 

The  normal  and  shear  stress  components  in  Cartesian  coordinates 
can  be  computed  from  the  equations. 


Gx  -  G^  cosV  +  %  s\no(  -t^  s\ri2<<  (ill) 


Sy    -   &n    sinV  H-  G*   coslc*  +  I|?  sin^oC  (15) 


T*y  =  1  (G|  -G7)  sin  2oi     -  T|J2  cos  <b<  (16) 


•where  cos  ©<f  - 


\  —  >  / 


sino(-   $[  s/n/7^  co5^  (l8) 

/cosh.^-  -cos?r? 

Equations  (l^),  (15)  and  (16)  are  derived  as  usual  from 
equilibrium  conditions  on  an  element  .(Fig.  lb) 
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BOUNDARY  CONDITIONS 


Equations  (11),  (12)  and  (15)  satisfy  the  required  conditions 
of  symmetry  and  antisymmetry  refered  to  in  the  previous  section. 
In  addition  to  this  G|  and  T^     must  vanish  at  ^-|o  .  Page  193,  i 

Therefore  at  each  point  on  the  elliptical  hole  we  must  satisfy 


&*>  =■   °    >     §  =  ?°  C1^) 


T!7-o    ;  t=to  (19b) 


For  the  case  of  one  elliptical  hole  in  an  elastic  finite 
flat  rectangular  plate  loaded  biaxially  we  must  satisfy  the  following 
two  conditions  (20)  and  (21 ),    in  addition  to  (19) 


Txx|  -   O  (20Q) 


6>j  -  G*     ;  .  V  _  ^y.  (20b) 


Tyy-  O  te,a) 


Gx  =  Gi      -   x  ^  ^L  (21b) 

2 
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Due  to  the  symmetry  about  both  the  x  and  y  axes  it  is 
enough  to  examine  only  one  quadrant  of  the  plate. 

Figure  2  illustrates  all  the  cases  under  consideration. 
The  numbers  appearing  in  parentheses  are  the  numbers  of  the 
equations  of  the  boundary  conditions  applied  for  each  particular 
case. 
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NUMERICAL  METHOD 

Initially  the  number  of  unknown  constants  appearing  in  equations 
(||)  up  to  (13)  was  determined  by  an  arbitrary  truncation  of  the  series. 

Later  the  behavior  of  the  series  was  examined  for  several  cases 
and  it  was  realized  that  the  series  doesn't  always  behave  well.  Similar 
convergence  problems  for  functions  involving  elliptical  coordinates  were 
pointed  out  in  reference  [6] 

In  this  case,  the  reason  is  believed  to  be  the  very  rapid  convergence 
of  the  coefficients  of  Dn  and  Hn  after  the  fifth  or  sixth  term,  resulting 
in  a  very  rapid  divergence  of  the  coefficients  of  Cn  and  Gn  for  all  of 
the  cases  that  were  examined. 

Finally  twenty  terms  plus  the  constant  term  of  the  series  were 
chosen. 

Using  the  Overdetermined  Collocation  method  (Ref.8  )  the  following 
cases  were  examined: 

a.-  Six  fixed  points  on  the  inner  boundary  and  another  six  fixed 
points  on  the  outer  boundary  were  selected  to  satisfy  the 
boundary  conditions.  (Fig. 2d) 
b.-  Decreasing  by  one  the  number  of  points  on  the  ellipse  with 
no  change  on  the  outer  edges  produced  no  difference  in  the 
numerical  results  obtained  in  case  (a), 
c-  With  six  points  on  the  ellipse  and  eight  on  the  outer 

boundary,  the  shape  of  the  curve  for  the  tangential  stress 
distribution  around  the  elliptical  hole  was  changed  to  a 
flatter  one  and  it  was  realized  that  the  boundary  conditions 


-  19  - 

were  not  satisfied  as  well  as  before  due  to  the  least 
squares  process  involved  in  the  numerical  method, 
d.-  The  case  of  five  points  on  the  inner  boundary  and  eight 
on  the  outer  one  was  examined  and  it  was  found  that 
the  error  introduced  was  smaller  than  in  case  (c),  but 
not  as  small  as  in  cases  (a)  and  (b).  (Verifying  the 
satisfaction  of  the  boundary  conditions) 
It  might  be  noted  here  that  a  smaller  number  of  points  than 
in  case  (a)  was  used  by  Christiansen  [  15  ]  for  a  similar  problem. 

For  the  case  of  axial  loading  parallel  to  the  minor  axis 
of  the  ellipse  and  for  values  of   3X<D.5  (wide  plate)  the 
expression  ^  approaches  infinity  and  no  results  could  be  obtained 
(Appendix  l). 

The  same  effect  happens  for  the  case  when  plates  with  3y>05 
are  loaded  parallel  to  the  major  axis  of  the  elliptical  hole. 

From  the  boundary  conditions  a  set  of  M  linear  algebraic 
equations  in  IT  unknowns  is  formed,  where  M>N  . 
In  matrix  notation 


AX  =  B  (22) 


where  A  is  the  M  by  N  matrix  of  the  coefficients  of 

the  unknown  constants  X. 

By  means  of  a  least-squares  method  the  normal  equations 
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are  obtained  and  the  system  in  matrix  notation  has  the  form 


T  T 


AAX  =  AB  (23) 


These  normal  equations  are  solved  by  elimination  using 
largest  pivotal  divisor  [Ref .  12  ] 

A  computer  program  was  developed  and  vrritten  in  FORTRAN  TV 
to  execute  all  the  required  steps  and  calculate  the  tangential 
stress  Sr^   around  the  ellipse  and  the  normal  stresses  along  the 
axes  of  symmetry  of  the  plate. 

Subroutines  for  matrix  operation  are  taken  from  reference  [12] 

Ail  the  necessary  calculations  for  the  problems  under 
consideration  are  carried  out  by  the  same  computer  program  by 
transfering  control  to  the  proper  boundary  conditions. 

The  required  input  information  for  each  plate  to  be 
examined,  must  be  typed  on  one  card  and  many  plates  can  be 
analyzed  the  same  time. 

Details  for  the  input  cards  are  given  on  comment  cards  in 
the  program. 

Additional  details  about  the  problem  formulation  are  given 
in  Appendix  1. 
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RESULTS 

Figures  3  up  to  31  are  typical  plots  of  results  obtained 
from  the  computer  program.  Whenever  possible,  experimental  data 
are  indicated.  Most  of  the  figures  are  self explanatory  but  additional 
information,  if  necessary,  is  given  in  the  discussion  of  each  figure. 

Figures  3  up  to  ]_1|  show  the  stress  distribution  around  the 
boundary  of  the  elliptical  hole  in  finite  plates  for  several  values 
of  ju.   and  1  . 

It  should  be  noted  that  in  cases  of  loading  parallel  to  the 
minor  axis,  the  tensile  stresses  predicted  are  2$  to  17$  higher 
than  experiments  reported  in  [  9]  ■while  the  much  smaller  compressive 
stresses  are  kCff0  to  60$  lower. 

When  the  external  loading  is  parallel  to  the  major  axis  of  the 
ellipse  the  agreement  with  [ 9 ]  lies  between  5$  and  25$. 

The  above  differences  in  compressive  stresses,  when  loading 
parallel  to  the  major  axis  are  of  minor  importance  from  an  engineering 
point  of  view  because,  these  stresses  are  very  small  compared  to  the 
maximum  tensile  ones. 

Figures  15  up  to  22  show  the  maximum  tensile  stresses  (SCF  and 
SCF')  as  function  of  1  for  several  values  of  yU.  .  (Finite  plate 
with  one  elliptical  hols). 

In  figure  23  several  cases  of  biaxial  loading  are  plotted. 
It  was  not  possible  to  examine  the  same  plate  used  in  [10]  because 
of  n     limitations  discussed  previously.  Other  plates  were  examined 
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and  the 'definitions  from  [10]  were  used. 

The  curves  a,ppear  to  have  the  same  shape. 

Figures  2*i  up  to  51  show  the  Gr>    and  (5/  distribution  around  the 
boundary  of  the  elliptical  hole  and  the  x  axis  respectively  for  typical 
uniaxial  and  biaxial  loading. 
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CONCLUSIONS 

The  stress  functions  assumed  have  been  shown  to  give  results 
that  are  comparable  with  available  experimental  data  and  therefore 
some  confidence  exists  in  the  theoretical  predictions  for  cases 
where  no  experimental  results  exist. 

The  stress  concentration  factors  for  plates  perforated  with 
elliptical  holes  under  in-plane  loading,  appear  to  be  much  higher 
in  comparison  with  similar  plates  perforated  with  circular  holes. 

Certain  disagreements  have  been  noticed  between  measured 
and  predicted  stresses  in  the  low  level  compressive  stress  range, 
while  the  predictions  for  the  high  level  tensile  stresses  are 
somewhat  higher  than  the  measured  ones  but  still  in  close  agreement. 
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RECOMMENDATIONS 

The  following  topics  are  recommended  for  further  investigation: 
a.-  Postulation  of  similar  stress  functions  in  terms  of 

hyperbolic  functions  instead  of  positive  exponentials,  in 
association  with  other  functions  decaying  rapidly  as 
increases  in  order  to  be  more  suitable  for  computer  use. 
Timpe  in  [13]  suggests  such  stress  functions, 
b.-i  Study  of  the  problem  for  cases  of  loading  that  produce 
ant i symmetrical  conditions  (shear  forces  applied  on  the 
edges  of  the  plate  for  example),  will  be  rather  easy  to  handle, 
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APPENDIX  1 


Equations  (ll),  (12)  and  (13)  may  be  rearranged  in  the  form  of  (2^), 
(25)  and  (26),  respectively.  Also  equations  (1*0,  (15)  and  (16)  are 
rearranged  in  the  form  of  (27),  (28)  and  (29)  respectively. 
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+  S~Gn[-L  (A s<n  -BsT,)sin^   -  ABsn  Cos^]     4- 

T7  =  2,4... 

oo 
+V"Hn[-L  (A4n-B4n)5in^  -Absn  COS  Tel   ]   4 

1l2,4... 
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In  the  above  expressions  for  £>*  ,  <5y  ;  T*y  : 
Ain   (i  =  1,  2,  3>  !0   represent  the  coefficients  of 

Cn,  Dn,  Gn,  Hn,  respectively, 
from  the  expression  for  Gj- 

Bin   (i  =  1,  2,  3,  it)   represent  the  coefficients  of 

Cn,  Dn,  Gn,  Hn,  respectively, 
from  the  expression  for  6y 

ABin  (i  =  1,  2,  3>  *0   represent  the  coefficients  of 

Cn,  Dn,  Gn,  Hn,  respectively, 
from  the  expression  for  Zgn 

The  value  of  _£"  is  obtained  by  using  the  formula 


i-=lioa  (j.±£ld-2\  (50 ) 

2   tU/-A'/«'   ' 
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It  is  obvious  that  as  1-rO     (i.e.  oc«W*)  b'-*a?  and  therefore 

^  — >  oo  , 

Analyzing  the  case  of  figure  2d,  for  example,  the  use  of 
equations  (2i[)  and  (26)  according  to  (19a)  and  (19b)  generates  a  set 
of  12  equations,  if  six  fixed  points  are  selected  on  the  boundary  of  the 
(quadrant)  ellipse 

Another  set  of  six  equations  is  generated  by  using  (27)  and  (29) 
according  to  (21),  if  three  points  are  selected  on  x=3. 

Finally  six  more  equations  are  generated  by  using  (28)  and  (29) 
according  to  (20),  if  another  three  points  are  selected  on  y=3. 

It  should  be  noted  that  the  number  of  unknown  constants  was 
selected  to  be  21. 

The  system  obtained  this  way  has  2)4  equations  with  21  unknowns 
and  it  is  solved  by  following  the  procedure  described  in  section 
NUMERICAL  METHOD. 

The  matrix  A^  is  formed  and  A  and  B  are  premulti plied  by 
AT  and  from  the  solution  of  the  system 


AT  A  X  =  AT  B 


we  obtain  the  values  of  the  unknown  constants. 

Summarizing  the  above  procedure,  we  have  chosen  a  number  of 
points  greater  than  the  number  of  unknown  constants  and  we  required 
that  the  boundary  conditions  be  satisfied  as  nearly  as  possible  at 
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these  points  (Least  square  sense)  Ref.  [  lH], 

Because  of  the  nature  of  the  assumed  series  solution,  the  number 
of  unknown  constants  is  always  odd  when  truncating  the  series  at  any 
value  of  n.  Therefore  the  system  of  equations  is  always  overdetermined 
at  least  by  one  equation. 
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